Angular Momentum Decomposition for an Electron 
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We calculate the orbital angular momentum of the 'quark' in the scalar diquark model as well as 
that of the electron in QED (to order a). We compare the orbital angular momentum obtained from 
the Jaffe-Manohar decomposition to that obtained from the Ji relation and estimate the importance 
of the vector potential in the definition of orbital angular momentum. 

PACS numbers: 



I. INTRODUCTION 

While the total angular momentum of an isolated system is uniquely defined, ambiguities arise when decomposing 
the total angular momentum of an interacting multi-constituent system into contributions from various constituents. 
Moreover, in a gauge theory, switching the gauge may result in shuffling angular momentum between matter and 
gauge degrees of freedom. In the context of nucleon structure, this gives rise to subtleties in defining these quantities 
that are more fundamental than those subtleties associated with the choice of factorization scheme. 

In the context of hadron structure, it is natural to perform a decomposition of the z component of the angular 
momentum as the I component of the quark spin has a partonic interpretation as a difference between parton densities. 
Indeed, in the light-cone framework, Jaffe and Manohar proposed a decomposition of the form [l| 
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whose terms are defined as matrix elements of the corresponding terms in the +12 component of the angular momentum 
tensor 

M+12 = ~ 2+75?+ + q + { f x ^) 2 q + + s + ' ij TrF +i A j + 2TrF +i (r x id) * A J . (2) 
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The first and third term in (|1I2[) are the 'intrinsic' contributions (no factor of fx) to the nucleon's angular momentum 
J z = +| and have a physical interpretation as quark and gluon spin respectively, while the second and fourth term 
can be identified with the quark/gluon orbital angular momentum (OAM). Here q + = ^7~7 + 9 is the dynamical 
component of the quark field operators, and light-cone gauge A + = A + A z = is implied. The residual gauge 
invariance is fixed by imposing anti-periodic boundary conditions A^(x^, oo~) = — Aj_(xj_, — oo~) on the transverse 
components of the vector potential. 

Since the quark spin term does not contain any derivatives, its manifest gauge invariance is evident. However, AG 
is also gauge invariant, as it is experimentally accessible. In gauges other than light-cone gauge, it is defined through 
a non-local operator [2]. The net parton OAM 

z z = Y. C «+ C l = \-\Y.^-\^ G ( 3 ) 
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can be related to differences between observables and is thus also obviously gauge invariant. However, similar to the 
case of AG, a manifestly gauge invariant operator defining C z would be non-local, reducing to a local expression in 
light-cone gauge only. For the individual OAMs the situation is more subtle and a detailed discussion can be found 
in Ref. 0. 

An alternative decomposition [3( of the nucleon spin 

\-\H^ + H LZ q + J Z g (4) 

9 9 
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FIG. 1: Schematic comparison between the two decompositions ([T]) and Q of the nucleon spin. In general, only |AS = i Aq 
is common to both decompositions. 



into quark spin, quark OAM, and gluon (total) angular momentum is obtained from the expectation value of 



M Oxy = J2 -q^ z q +^2qUfxiDYq 



fx E x B 



(5) 



with iD = id — gA. Its main advantages are that each term can be expressed as the expectation value of a manifestly 
gauge invariant local operator and that the quark total angular momentum J z = i Aq+L z can be related to generalized 
parton distributions (GPDs), using Q 



J z = - 



dx x [q{x) + E q (x, 0, 0)] 



(6) 



and can thus be measured in deeply virtual Compton scattering or calculated in lattice gauge theory. Its main 
disadvantage is that both quark OAM L z q as well as gluon angular momentum J z contain interactions through the 
vector potential in the gauge covariant derivative, which complicates their physical interpretation. 

Since the expectation value of q r y z 'S z q vanishes for a parity eigenstate, one can replace q^ z q — > q^ + Y, z q — q+jsq+j 
i.e. the Aq are common to both decompositions. This is not the case for all the other terms. For example, the angular 
momenta in these decompositions ((T]),(|][|) are not defined through matrix elements of the same operator and one should 
not expect them to have the same numerical value. However, no intuition exists as to how large that difference is. 

In the matrix element defining L z , one may make the replacement 



^ (f x iD^j q — g7 x iD^j q — > q (7 + Y~) (j x iDj q = q^OrX iD^j 



(7) 



provided that the expectation value is taken in a parity eigenstate. While the Dirac structure of the operator on the 
r.h.s. of ([7]) is now the same as that appearing in ([2]), Eq. ([7]) still contains the transverse component of the vector 
potential through the gauge covariant derivative, and therefore, even in light-cone gauge, C z and L z differ by the 

expectation value of q^_ (^f x gAj q + . While it has long been realized that in general C z ^ L z , The main purpose of 
this paper is to address this issue first in the context of a scalar diquark model and then in QED. 



II. ORBITAL ANGULAR MOMENTUM IN THE SCALAR DIQUARK MODEL 



In a two particle system we introduce center of momentum and relative _L coordinates as 



Pi 

r± 



= Pl_L d 

= x 1 r 1± 

= X2P1J 

= ri_L - 



P2J 



x 2 r 2 ± = xr 1± + (1 - x)r 2 ± 

- XlP2± = (1 - X)p%± - Xp 2 A 



r 2 j_ 



(8) 



(9) 
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where x\ — x and X2 = 1 — x are the momentum fractions carried by the active quark and the spectator respectively. 
For a state with Pi — 0, this implies pu_ = — P2L = kj_, allowing one to replace the OAM operator for particle 1 by 
(1 — x) times the relative OAM in such a state [J] 



L\ = ru. x pi ± = [R ± + (1 - x)tx] x kj 



(l-i)rixki = (l-i)r. (10) 
Likewise one finds that the 



Here we used that the internal wave function of a bound state satisfies (kj_) = 
expectation value of CJ^ can be replaced by the expectation value of xC z . 

We now use the above decompositions (PCJ) , (J3J) to calculate the OAM of the 'quark' in the scalar diquark model, 
where the two particle Fock space amplitudes read [5| 



(x,ki.) = 

T 2 

tp1 1 (x ) k±) = - 



M + - 

k 1 + ik 2 



>(*,ki) 
{x,k\) 



(11) 



with 



M 2 — - 



g/^l-x 

It' 2 4-TTl2 W 



Here g is the Yukawa coupling and M/m/X are the masses of the 'nu- 



cleon'/'quark'/diquark respectively. Furthermore x is the momentum fraction carried by the quark and k± = 
kj_ e — kj_ 7 represents the relative _L momentum. The upper wave function index j refers to the helicity of the 'nucleon' 
and the lower index to that of the quark. With the light-cone wave functions available (fTT|) . it is straightforward to 
compute either C z or J z , and hence IA from the Ji relation. 



This yields for the orbital angular momentum C z of the 'quark' 



L z = 



dx 



d 2 k ± 
16?r 3 



(1-x) 



Alternatively one may consider the OAM as obtained from GPDs using the Ji relation ([6]) as 

1 r 1 

Lg = - dx [xq(x) + xE(x, 0, 0) - Aq(x)} , 
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where 

xq{x) = Z8(\ — x 

Aq{x) = ZS(l-x 

xE(x,0,0) = 2Mg 2 x 



16tt 3 



1 '. 



d 2 k ± 
16tt 3 



(12) 



(13) 



(14) 



d 2 k. 



(1 - x) 2 (xm + M) 



Mg 2 



x(l - x) 2 {xm + M) 



16tt 3 [ x (\ - X )M 2 - (1 - x)m 2 - xX 2 - k 2 } 2 
As one may have expected, the wave function renormalization constant 



87r 2 -x{l - x)M 2 + (1 - x)m 2 + xX 2 ' 



Z = 1 



dx 



d 2 k ± 
16tt 3 







ij'l 



cancels in L z q , yielding 
LI = 



dx 



16tt 3 



(x-1) 



+ (ar+l) 



dxxE(x, 0,0). 



(15) 



(16) 



Since some of the above kj_ integrals diverge, a manifestly Lorentz invariant Pauli-Villars regularization (subtraction 
with heavy scalar A 2 — > A 2 ) is always understood. Evaluating the above integrals is tedious, but straightforward, and 
one finds 



CI = LI 



(17) 



as was expected since L z q in the scalar diquark model does not contain a gauge field term. However, there is no 
such identity for the OAM distribution. The distribution of the z component of the OAM C z (x) is defined as in 
(|12[) . but without the ar-integration. A comparison with (|13[) without x- integration, i.e. comparing C z (x) with 
L z (x) = \ [xq(x) + xE(x, 0, 0) — Aq(x)] (Fig. \2§ shows that, even in a model without gauge fields, L z (x) cannot be 
identified with the cc-distribution of C z for a longitudinally polarized nucleon Q • 
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FIG. 2: x distribution of the orbital angular momentum Cq(x) (full) compared to L z q (x) from the unintegrated Ji relation 
(dotted) in the scalar diquark model for parameters A 2 = 10m 2 = 10 A 2 . Both in units of 

III. ORBITAL ANGULAR MOMENTUM IN QED 

In QED, there are four polarization states in the Fock component. To lowest order, the respective Fock space 
amplitudes for a dressed electron with J 2 = +\ read 



*+i +1 (^k ± ) 
<J>Li +1 (z,k x ) 



x(l — a;) 
A; 1 + ik 2 
1-x 



^(x,ki) 



(18) 



(— - m) cA(a;,kj_) 



with J>(x,ki) = j(L 



M 2 - 



Using these light-cone wave functions, it is again straightforward to calculate the orbital angular momentum {T 
of the electron in the Jaffe-Manohar [l[ decomposition 



C ' : = 1 'l^ " 



*^ i _ 1 (x,k ± ) - 4 r L +1 (x,kx) 



(19) 



Likewise, it is straightforward to evaluate the OAM using the Ji relation 



1 f 1 

LZ e = o dx [xq e (x)+xE e (x,0,0)-Aq e {x)} 
* Jo 



(20) 
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with 0] 



xq e (x) 
Aq e (x) 
xE e (x,Q,0) = 4m 2 e 2 



Z5(l-x)+x 



ZS(l-x) + 
d 2 k 



d 2 k ± 
16tt 3 
2 k ± 



16tt 3 



x 2 (l-a;) 2 



(21) 



ir 2 (l - x) 2 



16tt 3 [ m 2(! 

Again the wave function renormalization constant 



X 2 x + k\f 



An 2 m 2 (l — x) 2 + X 2 x 



Z = 1 



dx 



d 2 k ± 
16tt 3 



drops out in (f^0|) . yielding 



r, z — - 

^"2 



(/.<: 



rf 2 k ± 

16tt 3 



(x-1) 



+ (x-l) 



+ (x + 1) 



V 1 i 



(22) 



1 f 1 

- / dra£; e (x,0,0). (23) 
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Because of the divergent integrals a Pauli-Villars subtraction with A 2 
end of the calculation, while A 2 3> m 2 . 

The evaluation of the above integrals is again straightforward, yielding 



A is understood and A 2 — ► at the 



and 



CI 



o 

Alt 



2tt 



- / ^(l-x 2 )log 



(1 ~x) 2 m 2 +xA 2 XZo 
(1 — x) 2 m 2 + xX 2 



a 

An 



A, A 2 2 



47T 



c£r(l + x 2 



4, A 2 



log 



(1 — x) 2 m 2 



(1 — x) 2 m 2 



(1 - ir) 2 m 2 + xA 2 
(1 — x) 2 m 2 + a;A 2 (1 — x) 2 m 2 + xA 2 (1 — x) 2 m 2 + xA 2 



(24) 



(25) 



Both C z and L z e are negative, regardless of the value of A 2 (as long as A 2 > A 2 ). In the case of C z e the physical reason 
is helicity retention [6j, which favors the emission of photons with the spin parallel (as compared to anti-parallel) 
to the original quark spin — particularly when x — > — resulting more likely in a state with negative OAM. The 
divergent parts of C\ and L\ are the same so that their difference is UV finite (Fig. [3]) 



A^oo 



Att 



CI -LI 

Applying these results to a (massive) quark with J z = +i yields to 0(a s ) 

CI 



(26) 



(27) 



i.e., for a s as 0.5 about 10% of the spin budget for this quark. 

In QCD, the gluon spin is experimentally accessible, but the gluon OAM C z is not. On the other hand, the gluon 
(total) angular momentum J z appearing in the Ji decomposition is accessible, either indirectly (by subtraction, using 
quark GPDs from lattice QCD and/or DVCS), or directly, using by calculating gluon GPDs on a lattice and/or deeply 
virtual J/ip production. Even though \AG and J z belong to two incommensurable decompositions of the nucleon 
spin, one may thus be tempted to consider the difference between these two quantities, hoping to learn something 
about gluon OAM. Subtracting |T]) from Q, it is straightforward to convince oneself that 



Jn 



-^G = C z g 



(28) 
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FIG. 3: Cutoff dependence of L% (full) and L% (dotted). Both in units of 



i.e. numerically J| — \ AG differs from C z by the same amount that C z q differs from ^ L q . In our QED example, 
with 



A7 



dx 



16tt 3 



being the photon spin contribution, one thus finds (for A — > 0, A — » 00) 

4-^ = ^ + £- 



(29) 



(30) 



As was the case in P5|) . ^jp appears to be a small correction, but one needs to keep in mind that for an electron J z , 
A7, and C z are also only of order a. 



IV. DISCUSSION AND SUMMARY 



We have studied both the Jaffe/Manohar, as well as the Ji decomposition of angular momentum in the scalar 
diquark model, as well as for an electron in QED to order a. As expected, both decompositions yield the same 
numerical value for the fermion OAM in the scalar diquark model, but not in QED. This calculation demonstrates 
explicitly that the presence of the vector potential in the manifestly gauge invariant local operator for the OAM does 
indeed contribute significantly to the numerical value of the OAM. While the numerical value for difference between 
the fermion OAM in these two decompositions in QED appears to be small (■¥-), one should keep in mind that the 
OAM itself is of the same order a. Moreover, applying the same calculation to a massive quark in QCD yields a 
contribution from the vector potential term to the angular momentum of the quark of about —10% (for as ~ 0.5). 

The sign of the contribution to the angular momentum arising from the vector potential is also significant in light 
of recent lattice results for the contributions from the u and d quark OAM to the nucleon spin [|[, yielding L z u < 
and > 0. The signs of the lattice results are thus exactly opposite to what one would have expected on the basis 
of relativistic quark models, such as the bag model, where the OAM arises from the lower Dirac component and its 
expectation value is thus positively correlated to the expectation value of the quark spin. While the lattice results 
still neglect insertions of the operator into disconnected quark loops, this does not affect L z u — L%, and the sign of 
that difference should be reliable. In Ref. @, evolution has been proposed to explain this apparent discrepancy, as a 
quark acquires OAM in the direction opposite to its spin from virtual gluon emission (see Fig. [3]) . Our result adds to 
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that effect in the sense that the vector potential also adds a contribution to the OAM that is in the opposite direction 
from the quark spin. Such a shift would imply L z u > L z u and C z d < L z d , moving C z q closer to the quark-model-based 
intuitive expectation than L z . 
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